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' By means of the adiabatic self-consistent collective coordinate method and the pairing 

plus quadrupole interaction, we have obtained the self-consistent collective path connecting 
^s^j , the oblate and prolate local minima in 68 Se for the first time. Result of calculation indi- 

■ cates importance of triaxial deformation dynamics in the oblate-prolate shape coexistence 

' phenomena. 

m ' 
o 

l/-) ■ Shape coexistence phenomena are typical examples of large amplitude collective 

motion in nuclei. These phenomena imply that different solutions of the Hartree- 
Fock-Bogoliubov (HFB) equations (local minima in the deformation energy surface) 
appear in the same energy region and that the nucleus exhibits large amplitude collec- 
tive motion connecting these different equilibrium points. The identities and mixings 
of these different shapes are determined by the dynamics of such collective motion. 
Some years ago, we have proposed a new method of describing such large-amplitude 
collective motion, which is called Adiabatic Self-Consistent Collective Coordinate 
(ASCC) methodP It yields a new method of solving the basic equations of the SCC 
methocP* using an expansion in terms of the collective momentum. It does not as- 
sume a single local minimum, so that it is expected to be suitable for the description 
of the shape coexistence phenomena. The ASCC method also enables us to include 
the pairing correlations self-consistently, removing the spurious number fluctuation 
modes. To examine the feasibility of the ASCC method, we have first applied it 
to an exactly solvable model called the multi-0(4) model, which is a simplified ver- 
sion of the pairing-plus-quadrupole (P+Q) interaction modelP It is shown that the 
method yields a faithful description of tunneling motion through a barrier between 
the prolate and oblate local minima in the collective potential. 4 ' 

In this Letter, we give a brief report of our first application of the ASCC method 
to a realistic P+Q interaction model. We illustrate its practicality taking as a typical 
example the oblate-prolate shape coexistence phenomenon in 68 Se recently observed 
in experiments.^ The self-consistent collective path obtained successfully by means 
of the ASCC method is found to run approximately along the valley connecting 
the oblate and prolate local minima in the collective potential energy landscape. 
To the best of our knowledge, this is the first time that a self-consistent collective 
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path is obtained for realistic situation starting from the microscopic P+Q Hamil- 
tonian. We note that a similar approach to large amplitude collective motions was 
recently pursued by Almehed and Waletj^ although they discussed different nuclei 
and encountered some difficulties in obtaining self-consistent collective paths. 

We assume that large-amplitude collective motions are described by a set of 
time-dependent HFB state vectors \cj)(q,p,ip,N)} parametrized by a single collective 
coordinate q, the collective momentum p conjugate to q, the particle number N and 
the gauge angle <p conjugate to N. As discussed in Ref. P), the state vector can be 
written 

\<Kq,p,<p,N)) =e-*e*\<f>(q,p,N)) =e^^\<l>(q)). (1) 

Making an expansion with respect to p and requiring that the time-dependent vari- 
ational principle be fulfilled up to the second order in p, we obtain the following set 
of equations to determine |0(g)), the infinitesimal generator Q(q), and its canonical 
conjugate P(q): 

HHq)\H M (q)\4>(q))=0, (2) 

5 (0(g)| [H M (q),Q(q)] ~ ~B(q)P(q) |0(g)) = 0, (3) 

i 

5 (0(g) | [H M (q), lp{q)]-C{q)Q{q)-^--[{H M {q), (H-X(q)N) A ],Q(q)]-^N |0(g)) = 0. 

(4) 

Here 

H M (q) = H-X(q)N-^Q(q) (5) 
is the Hamiltonian in the moving frame; 

_ d 2 V 1 dBdV 

C{q) -W + Wq)^q^q {) 

is the local stiffness; (H — XN)a represents the two-quasiparticle creation and anni- 
hilation parts of (H — XN); Q(q) and P(q) satisfy the canonical variable condition 

(Hq)\[Q(q),P(q)]\<P(q))=i- (7) 

Once |0(g)) and the infinitesimal generators are determined for every values of g, we 
obtain the collective Hamiltonian TC(q,p) = \B(q)p 2 +V {q) with the collective poten- 
tial V(q) = (0(g) | H |0(g)) and the inverse mass B(q) = - (0(g) | [{H, Q{q)],Q(q)) |0(g)). 

Table I. Spherical single-particle orbits and their energies used in the calculation. 
Energies relative to those of lgg/2 are written in MeV. 



orbits 


I/7/2 


2P3/2 


I/5/2 


2fl/2 


I39/2 


2^5/2 


lff7/2 


3Si/ 2 


2d 3 /2 


protons 


-8.77 


-4.23 


-2.41 


-1.50 


0.0 


6.55 


5.90 


10.10 


9.83 


neutrons 


-9.02 


-4.93 


-2.66 


-2.21 


0.0 


5.27 


6.36 


8.34 


8.80 
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We use the P+Q interaction model with the prescriptions of Ref. IB} for the 
microscopic Hamiltonian H, but here the pairing and quadrupole force parameters 
are chosen as G = 0.320 MeV (for both protons and neutrons) and x' = 0.248 
MeV so that the constrained HFB potential energy surface (shown by contour lines 
in Fig.l) exhibits two local minima at prolate and oblate shapes, whose pairing 
gaps, quadrupole deformation and energy difference approximately reproduce those 
obtained in a recent Skyrme-HFB calculation by Yamagami et al. 7) The spheri- 
cal single-particle energies are taken from those of the modified oscillator model of 
Ref.|SJ) and listed in Table I. In this way the effective Hamiltonian provides a suitable 
situation with which shape coexistence dynmamics can be studied, although further 
improvements, e.g., by including the quadrupole pairing and/or neutron-proton pair- 
ing interactions, may better be taken into account for quantitative comparison with 
experimental data. 

We have used the following algorithm to solve the set of ASCC equations ©, 
(jHJ), @ and (0). Let the state vector \4>(q)} be known at a specific value of q. We first 
solve the local harmonic equations in the moving frame (the moving frame RPA), 
|2J) and (JIJ, under the condition (J7J) to obtain Q(q) and P{q). This is done by a 
straightforward extension of the procedure described in Ref. QJ. We then construct 
a state vector at the neighboring point q + 5q by using the infinitesimal generator 
P(q) as 

\ttq + Sq)) = er i5 * p W\<l>(q)), (8) 

and solve the moving frame RPA with respect to this state to obtain Q{q + 5q) 
and P(q + 5q). Though the above \(p(q + 5q)) does not necessarily satisfy the HFB 
equation in the moving frame ©, we can use this state vector as an initial solution 
of J2J) at q + 5q. We search for the solution of © under the constraints 

{<l>(q + 8q)\N\<t>(q + 6q)) = N, (9) 
(<t>(q + 6q)\Q(q)\<t>(q + 6q))=6q (10) 

by means of the gradient method. Here the nucleon-number constraint © is actually 
applied for both proton and neutron numbers. Equation ()10() is the constraint for 
the increment 5q of the collective coordinate. After finding a solution of Eq. @, we 
renew Q(q + Sq) and P(q + 5q) by solving again the moving frame RPA equations, © 
and @, for the new state vector \(j>(q + 5q)) obtained above. Then we again solve Eq. 
© with the renewed Q{q + 5q). If the above iterative procedure converges, we get 
the self consistent solutions that satisfy Eqs. ©, ©, @ and (J7J) simultaneously at 
q + Sq, and we can proceed to the next point q + 25q. In actual numerical calculation, 
we start the procedure from one of the HFB local minimum and examine whether 
we arrive at the other local minimum by going along the collective path obtained 
above. We have checked that the same collective path is obtained by starting from 
the other local minimum and proceeding in an inverse way. 

Carrying out the above procedure we have successfully obtained the collective 
path connecting the oblate and prolate local minima in 68 Se. The result is shown in 
Fig. 1. The deformation parameters (3 and 7 are here defined as usual through the 
expectation values of the quadrupole operators^ Roughly speaking, the collective 
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Fig. 1. The bold curve represents the ASCC 
path connecting the oblate and prolate 
minima in 68 Se projected on the (/3, 7) 
plane. The contour lines are calculated by 
the conventional constrained HFB method 
and plotted for every 100 keV. 
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Fig. 2. Collective potential V(q) plotted as 
a function of the collective coordinate q. 
Here the origin of q is chosen to coin- 
cide with the prolate local minimum and 
its scale is defined such that the collective 
mass M(q) = 1 MeV -1 . 
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Fig. 3. Collective Mass M(s) with respect to 

the geometrical length s along the collec- Fig. 4. The triaxiality parameter 7 plotted as 
tive path in the (/3, 7) plane is plotted as a a function of the collective coordinate q. 

function of the triaxiality parameter 7. 



path goes through the valley that exists in the 7 direction and connects the oblate 
and prolate minima. If one treats the (3 as collective coordinate and connects the 
oblate and prolate shapes through the spherical point, variation of the potential 
energy would be much greater than that along the collective path we obtained. The 
potential energy curve V(q) along the collective path evaluated by the ASCC method 
is shown in Fig. 2. We have defined the scale of the collective coordinate q such that 
the collective mass M(q) = B(q)- 1 = 1 MeV -1 . The collective mass as a function of 
the geometrical length s along the collective path in the (/?, 7) plane may be defined 
by M(s) = M(q){ds/dq)~ 2 with ds 2 = d(3 2 + f3 2 d-y 2 . This quantity is presented 
in Fig. 3 as a function of 7. The triaxial deformation parameter 7 is plotted as a 
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Fig. 5. Neutron and proton pairing gaps, A n 
and A p , plotted as functions of 7. 
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6. Lowest three eigen- frequencies squared, 
ui 2 — BC, of the moving frame RPA, plot- 
ted as functions of 7. 



function of q in Fig. 4. Variations of the pairing gaps and of the lowest few eigen- 
frequencies of the moving frame RPA along the collective path are shown in Figs. 5 
and 6. The solid curve in Fig. 6 represents the squared frequecy u> 2 (q) = B(q)C(q), 
given by the product of the inverse mass B(q) and the local stiffness C(q), for the 
moving frame RPA mode that develops from the 7-vibration in the oblate and prolate 
limits and determines the infinitesimal generators Q(q) and P(q) along the collective 
path. The other two curves are solutions of the moving frame RPA having characters 
of the collective rotational motion and the (3- vibration, which are however irrelevant 
to the collective path. Note that the frequency of the 7 mode becomes imaginary in 
the region 12° < 7 < 45°. These results will reveal interesting dynamical properties 
of the shape coexistence phenomena in 68 Se. For instance, the large collective mass 
in the vicinity of 7 = 60° (Fig. 3) might increase stability of the oblate shape in the 
ground state. Detailed discussions on these quantities as well as the solutions of the 
collective Schrodinger equation will be given in a forthcoming full-length paperP) 

In summary, we have applied the ASCC method to the oblate-prolate shape co- 
existence phenomena in 68 Se. It was found that the collective path goes through the 
valley of the potential energy landscape in the (/?, 7) plane along which the triaxial 
deformation parameter 7 changes between 0° and 60° keeping the axially symmet- 
ric deformation parameter (3 approximately constant. This is the first time that a 
self-consistent collective path between the oblate and prolate minima is obtained for 
the realistic P+Q interaction model. Currently, the generater coordinate method 
has often been used to describe variety of shape coexistence phenomena taking the 
(3 as the generater coordinated The triaxial shape vibrational degrees of freedom is 
ignored also in the extensive microscopic calculation of Ref . I1U|) . The result of the 
ASCC calculation, however, strongly indicates the necessity of taking into account 
the 7 degree of freedom at least for describing the oblate-prolate shape coexistence 
in 68 Se. Effects of triaxial deformation dynamics on various properties of shape co- 
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existence, including results of calculation for neighboring nuclei, will be discussed in 
a full-length paperPJ 
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